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\  Abstract 

~^This  investigation  was  undertaken  to  determine  the 
sensitivity  of  the  damping  factor  in  a  large  space  structure 
CLSS)  to  small  changes  in  nonstructural  mass  and  structural 
element  stiffness.  Revision  3  of  the  ACOSS  2  model,  devel¬ 
oped  by  the  Charles  Stark  Drapes  Laboratory,  Inc.  (CSDL) , 
was  used  as  the  model  of  the  LSS.  Various  combinations  of 
the  mirror  masses  in  this  large  space  telescope  were  varies 
by  up  to  lQ7o,  selected  structural  elements  were  stiffened 
by  increasing  their  cross-sectional  areas  by  10  and  50%, 
and,  finally,  two  structural  elements  in  the  middle  of  the 
telescope  were  stiffened  to  represent  the  addition  of  a 
lumped  mass  located  away  from  the  control  system  sensors 
and  actuators.  A  control  system  of  21  collocated  sensors 
and  actuators,  positioned  at  the  top  and  bottom  of  the 
telescope,  was  used  in  this  analysis. 

The  analysis  was  accomplished  using  NASTRAN  for  the 
finite  element  analysis,  and,  after  selecting  certain 
vibration  modes  for  further  study,  the  conq>lex  conjugate 
pairs  used  to  determine  the  damping  factors  were  calculated. J: 
This  last  step  was  accomplished  through  spillover  reduction 
and  reduction  of  the  closed  loop  matrix  to  a  lower  tri¬ 
angular  form.  This  procedure  was  repeated  for  each  of  the 
46  perturbations  applied  to  the  model  in  this  research. 


The  results  indicate  that  the  damping  factors  for  the  CSDL 
model  are  quite  sensitive  to  small  perturbations  and  that 
it  is  very  difficult  to  predict  the  effect  that  a  pertur¬ 
bation  to  the  model  will  have  on  the  damping  factors.  If 
the  vibrations  in  an  LSS  are  to  be  altered  by  changing^masT^ 
or  stiffness ,  a  thorough  analysis  is  required  to  ensure  that 
increasing  the  damping  for  one  mode  does  not  result  in  an 
undesirable  change  in  the  damping  for  another  mode . 
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THE  EFFECT  OF  MASS  AND  STIFFNESS  CHANGES  ON  THE 


DAMPING  FACTOR  IN  A  LARGE  SPACE  STRUCTURE 
AS  REPRESENTED  BY  THE  CSDL  2  MODEL 


I,  Introduction 


Back^roimd 


'^Deployment  of  large  space  structures  for  commu¬ 
nications.  space  defense,  power  generation,  manufac¬ 
turing,  and  research  has  become  a  major  objective  of 
space  planners  for  the  mid-1980's  and  beyond.  These 
systems  typically  combine  large  size  with  extremely 
rigorous  pointing  and  surface  figure  performance 
requirements."  (Ref  1:1-1) 

In  addition  to  these  performance  requirements ,  the  mass  of 
a  large  space  structure  (LSS)  must  be  kept  as  small  as  pos¬ 
sible  to  minimize  launch  costs  and  the  structure  must  be 
designed  to  be  easily  assembled  once  delivered  to  the  proper 
orbit.  These  constraints  generally  "result  in  highly  flex¬ 
ible  spacecraft  which  exhibit  poor  dimensional  precision. 

Our  current  ability  to  conceptualize  space  systems,  to  se¬ 
lect  structural  materials,  and  to  build  mechanisms  for  on- 
orbit  construction  and  deployment  has  now  far  outstripped 
our  demonstrated  ability  to  control  the  resulting  structure." 
(Ref  1;1-1) 


In  order  to  study  the  structural  characteristics  of  an 
LSS,  the  Charles  Stark  Draper  Laboratory,  Inc,  (CSDL)  of 
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Cambridge,  MA  developed  a  model  representative  of  an  above 
described  structure.  This  model,  known  as  CSDL  2,  is  de¬ 
picted  in  Figure  1.  It  is  a  large  space  telescope  measur¬ 
ing  28  meters  tall  and  20  meters  wide,.  The  beam/ truss 
structure  is  made  of  epoxy -graphite  tubes  varying  in  diam¬ 
eter  from  20  to  40  cm.  Including  the  mirrors,  solar  panels, 
and  equipment  package,  the  telescope  weighs  about  9300  kg. 
Three  desired  features  were  incorporated  in  the  design  of 
this  model;  1)  a  structural  design  based  on  reasonable 
sizes  and  weights,  2)  a  simple,  vmclassified  optical  system 
with  associated  performance  measures  and  weights,  and  3)  a 
set  of  disturbances  typical  of  equipment  vibration  and  at¬ 
titude  control,  (Ref  2:1)  This  model  is  being  analyzed  in 
the  Active  Control  of  Space  Structures  (ACOSS)  program, 
sponsored  by  the  Defense  Advanced  Research  Projects  Agency 
(DARPA) ,  to  determine  the  adequacy  of  various  active  struc¬ 
tural  control  methods  xmder  development.  Another  effort, 
the  Vibration  Control  of  Space  Structures  (VCOSS)  program, 

"is  intended  to  study  the  application  of  ACOSS  technology 
to  an  actual  spacecraft  design",  to  include  specifying  the 
actual  control  hardware  to  be  used  and  assessing  the  control¬ 
lability  of  the  structure  using  this  hardware.  The  minimum 
mass  VCOSS  model  is  based  on  the  ACOSS  model  but  has  been 
modified  to  '^reduce  the  structural  mass  to  the  minimum  re¬ 
quired  to  maintain  structural  stability,  This  system  will 
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CSDL  2  Conceptual  Model 


rely  entirely  on  the  control  system  for  vibration  supression 
to  meet  the  performance  requirements."  (Ref  3:1)  If  hard¬ 
ware  requirements  are  to  be  specified,  the  system  model  must 
accurately  represent  the  vibrational  properties  of  the  actual 
structure.  In  this  reduced  mass  structure,  the  Interactions 
between  the  mirrors  and  flexible  support  structure  and  be¬ 
tween  the  equipment  section  and  optical  support  truss  are 
the  most  important  properties.  Thus,  the  modeling  of  these 
areas  was  given  particular  attention  in  the  VCOSS  program, 
resulting  in  revision  three  of  the  CSDL  2  model.  This  is 
the  model  that  is  analyzed  in  this  thesis.  The  wall  thick¬ 
ness  of  the  tubes  in  revision  three  varies  from  ,03  to  ,067 
cm,  the  radius  ranges  from  3,6  to  8,1  cm,  and  the  total 
weight  is  down  to  89.63  kg.  This  model  is  thoroughly  dis¬ 
cussed  in  reference  3. 

The  ACOSS  FIVE  research,  study  has  indicated  that  con¬ 
trol  system  performance  associated  with  detailed  structural 
models  "can  be  extremely  sensitive  to  changes  in  the  struc¬ 
tural  model"  and  these  sensitivities  may  result  from  the 
"detailed  way  a  portion  or  all  of  the  structure  is  built," 
(Ref  1:4-1), 

These  research  programs  use  the  NASTRAN  finite  element 
structural  analysis  computer  program  developed  by  NASA  and 
released  to  the  public  in  1969,  The  NASTRAN  program  re¬ 
quires,  as  part  of  its  input  data,  stiffness  and  inertia 
properties  for  each  structural  element.  Non -structural 


4 


masses,  such  as  the  mirrors  on  the  CSDL  2  model,  are  handled 
as  point  masses,  that  is,  their  mass  is  assumed  to  be  con¬ 
centrated  at  specific  points  on  the  structure.  These  mass 
distributions  are  also  inputs  to  the  NASTRAN  program. 

NASTRAN  can  provide,  among  other  things,  data  on  the  natural 
vibration  frequency  and  mode  shape  associated  with  each 
structural  element.  CRef  4»316)  The  natural  frequency  is 
the  frequency  at  which  the  structure  will  vibrate  if  there 
are  no  external  forces  being  applied  to  it.  The  shape  that 
each  structural  element  assumes  when  it  is  vibrating  is  the 
mode  shape  of  the  element.  (Ref  SilQO).  Each  structural 
element  will  have  specific  natural  frequencies  and  mode 
shapes  associated  with  it. 

Control  system  design  starts  with  the  finite  element 
structural  model,  which  yields  natural  frequencies  and 
mode  shapes.  The  control  system  engineer  then  selects  ap¬ 
propriate  actxiators ,  which  produce  forces  or  torques ,  and 
sensors,  used  to  measure  the  amount  and/or  rate  of  displace¬ 
ment,  The  next  step  in  the  control  design  process,  the 
determination  of  where  to  locate  the  sensors  and  actuators , 
was  addressed  by  the  ACOSS  FIVE  report,  (Ref  1}4-15)  The 
desired  result  of  control  design  is  a  control  system  of 
sensors  and  actuators  that  is  able  to  effectively  negate 
an  induced  vibration  without  exciting  a  vibration  in  the 
rest  of  the  structure.  The  set  of  natural  frequencies  and 
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mode  shapes  associated  with  a  structure  (hereafter  called 
the  ({>  matrix)  can  be  analyzed  in  conjunction  with  the  de¬ 
sired  location  and  orientation  of  the  sensors  and  actuators 
(called  the  D  matrix)  to  determine  the  effectiveness  of  the 
control  system.  The  product  of  the  transpose  of  the  ^  ma¬ 
trix  and  the  D  matrix  gives  the  B  matrix,  whose  rows  phys¬ 
ically  correspond  to  the  amplitude  of  each  mode  at  the 
various  actuator  locations.  If  the  sensors  and  actuators 
are  collocated,  then  the  C  matrix,  x<7hose  columns  physically 
correspond  to  the  amplitude  of  each  mode  at  the  various 
sensor  locations,  is  the  transpose  of  the  B  matrix. 

(Ref  6:2)  The  B  and  C  matrices  can  be  analyzed  to  deter¬ 
mine  the  damping  factor  of  the  structure,  which  represents 
a  measure  of  the  total  amount  of  damping  (both  active  and 
passive)  present  in  the  system.  CRef  7:62)  There  will  be 
a  certain  amount  of  damping  in  the  structure  that  comes 
from  its  resistance  to  bending  and  twisting.  The  control 
system  is  designed  to  increase  the  damping  factor  and  hence 
cause  the  magnitude  of  the  vibration  to  decrease  faster  than 
it  would  in  the  absence  of  any  active  control. 

The  development  of  the  initial  CSDL  2  model  was  com¬ 
pleted  by  the  Draper  Laboratory  in  May  of  1980.  The  product 
of  this  effort,  the  space  telescope  model,  is  an  example  of 
a  typical  large  space  structure  that  can  be  used  to  develop 
and  evaluate  control  systems  for  1)  positioning  the  struc- 
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ture,  and  2)  controlling  structural  vibrations  caused  by 
thrustors  (actuators)  and  onboard  equipment.  As  discussed 
in  section  five  of  the  ACOSS  Model  2  research  report  (Ref  2) 
two  sets  of  parameter  variations  on  the  CSDL  2  model  were 
input  to  the  NASTRAN  program  to  '^assess  the  sensitivity  of 
the  control  system  to  changes  in  the  natural  frequencies 
and  mode  shapes,'*  (Ref  2; 45)  The  selection  of  these  vari¬ 
ations  in  stiffness  and  mass  distribution  represent  "two 
point  designs  from  an  infinite  nijmber  of  possible  variations 
and  "In  order  to  fully  evaluate  the  stability  limits  of  a 
particular  control  system  design,  it  may  be  necessary  to 
evaluate  several  perturbed  models."  (Ref  2:45) 


Objectives  of  the  Research 

The  objective  of  this  research  is  to  expand  the  vibra¬ 
tional  analysis  of  a  typical  LSS,  represented  by  the  CSDL  2 
model,  with  respect  to  variations  in  beam  stiffness  and  non- 
structural  mass.  This  research  will  determine:  1)  the 
effect  that  these  variations  have  on  the  controllability 
of  the  structure  by  looking  at  the  change  in  the  damping 
factor  for  selected  modes  of  vibration,  2)  the  particular 
structural  elements  or  non-structural  masses  that,  when 
vAried  slightly,  have  a  significant  effect  on  the  magnitude 
of  the  damping  factor,  and  3}  the  feasibility  of  adding 
non-^structural  mass  or  stiffening  structural  elements  to 
increase  the  controllability  of  an  existing  LSS.  This 


informaCion  can  be  used  by  structural  design  engineers  to 
specify  allowable  manufacturing  tolerances  for  certain  struc¬ 
tural  elements  of  an  LSS  and  by  control  engineers  to  develop 
or  improve  control  systems  capable  of  handling  the  wide 
range  of  vibrations  in  an  LSS, 

Methodology  Overview 

The  analysis  of  the  effect  of  small  variations  in 
stiffness  and  mass  is  accontplished  using  five  computer  pro¬ 
grams.  Selected  perturbations  are  applied  to  the  CSDL  2 
model  and  then  NASTRAN  is  run  to  determine  the  first  30 
frequencies  and  mode  shapes .  Program  FORMAT  reformats  the 
NASTRAN  output  for  input  to  program  SELECT,  SELECT  picks 
off  the  particular  modes  to  be  analyzed  C12  modes  in  this 
research)  and  formats  the  data  for  input  to  ZETA,  In  pro¬ 
gram  ZETA,  the  B  and  C  matrices  are  calculated  from  the  ((> 
matrix  and  the  D  matrix,  and  the  output  is  formatted  for 
input  to  DOFB34,  D0FB34  CDirect  Output  Feedback  34) , 
through  spillover  reduction,  reduces  the  closed  loop  equa¬ 
tion  matrix  to  a  lower  triangular  form  and  produces ,  among 
other  things,  two  eigenvector  pairs  per  mode.  C^ef  8:49) 
These  pairs  are  used  to  calculate  the  damping  factor  for 
the  mode;  the  damping  factor  is  then  used  as  the  performance 
index  to  determine  if  the  structure  is  becoming  more  or  less 
stable  as  a  result  of  variations  in  stiffness  and  mass. 


A  detailed  discussion  on  the  use  of  these  programs  is 


provided  in  Appendix  B.  The  equations  used  in  ZETA  and 
the  calculating  of  the  damping  factor  are  discussed  in 
Chapter  3.  DOFB34  is  developed  and  analyzed  in  Reference 
8  and  will  not  be  extensively  discussed  in  this  report. 
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II.  The  CSDL  2  Model 


Basic  Model 

The  CSDL  2  Model  CFigure  1)  represents  a  three -mirror , 
wide-angle  optical  space  telescope.  The  two  major  compo¬ 
nents  of  this  structure  are  the  optical  support  structure 
and  the  equipment  section,  connected  to  each  other  by  springs 
at  three  points  to  allow  either  active  or  passive  vibration 
isolation.  The  optical  support  structure,  consisting  of 
the  upper  mirror  support  truss,  the  lower  mirror  support 
truss,  and  the  metering  truss,  contains  the  four  optical 
surfaces.  The  upper  mirror  support  truss  contains  the 
primary  and  tertiary  mirrors,*  the  secondary  mirror  and  the 
focal  plane  are  contained  in  the  lower  mirror  support  truss, 
The  metering  truss  maintains  mirror  separation  and  is  the 
key  section  when  examining  the  focus  of  the  telescope,  The 
equipment  section  is  modeled  as  a  central  rigid  body  with 
two  flexible  solar  panels  cantilevered  from  it,  and  is 
assumed  to  contain  all  of  the  guidance,  navigation,  control, 
and  power  systems  required  for  the  operating  of  the  tele¬ 
scope.  The  full  structure  is  approximately  28  meters  high 
and  weighs  8963  kg  in  the  VCOSS  strength^controlled  design 
that  is  analyzed  in  this  research. 

The  finite  element  model  of  the  structure  CFigure  2) 
contains  59  node  points,  31  of  which  represent  nodes  where 
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structural  elements  connect.  The  other  8  nodes  were  added 
to  more  accurately  model  the  mirrors  and  equipment  section 
(Ref  3) .  The  coordinates  of  each  node  are  contained  in  the 
NASTRAN  input  deck  CAppendix  A)  under  the  title  of  GRID 
points.  The  location  of  the  node  points  is  shown  in  Figure 
2. 

There  are  137  structural  elements  in  the  CSDL  2  model 
made  of  epoxy<-graphlte  tubes  varying  from  3.6  to  8,1  cm  in 
radius.  The  CBAR  entries  in  the  NASTRAN  input  deck  list 
the  node  points  connected  by  each  structural  element.  The 
optical  support  structure  is  designed  to  act  as  a  truss , 
but  it  is  assumed  that  all  joints  allow  for  a  full  moment 
connection.  Thus,  both  bending  and  axial  stiffness  are  in> 
eluded  for  all  elements.  (Ref  2:2)  The  mass  of  each 
structural  element  is  lumped  at  its  node  points  and  is 
computed  by  NASTRAN  using  the  cross-sectional  area  and 
length  of  each  member  and  a  material  density  of  1720  kg/m^. 
For  this  VCOSS  model,  the  minimum  allowable  wall  thickness 
for  the  epoxy->graphite  tubes  is  .03  cm,  and  the  ratio  of 
radius  to  thickness  (r/t)  was  defined  to  be  a  maximum  of 
120.  The  length  to  radius  ratio  (l/r)  is  defined  to  be  a 
maximum  of  400.  The  natural  frequency  of  each  member  "was 
constrained  to  be  greater  than  10  hz  to  prevent  significant 
interaction  between  systems  vibrations  and  local  vibrations." 
(Ref  3jl7)  Thus,  the  radius  of  each  tube  is  the  maximum 


12 


value  of  the  three  radius  values  determined  by  the  above 
criteria.  (Ref  9)  That  is,  the  radius  r  is  the  maximum  of; 


1) 

r  = 

r/t  (tmin)  «  120  (.0003m)  =  .036  meters 

(1) 

2) 

r  = 

1  ®  (.0025)1  meters 

(2) 

3) 

r  * 

1,1104  X  10"^^  tiil^  meters 

(.3) 

where  oi  ■»  2TrClQ  hz)  radians/ second  and  1.1104E»“05  represents 
the  interaction  between  the  local  buckling  load,  1/r,  Youngs 
modulus  E,  area,  and  a  constant,  K,  that  accounts  for  the 
end  conditions.  Additionally,  r  is  rounded  to  the  nearest 
.0002  m. 

The  cross"«'Sectional  area  and  moments  of  inertia  for 
each  of  the  structural  elements  are  listed  in  the  NASTRAN 
data  deck  under  the  title  PEAR.  In  developing  this  model, 
the  following  equations  were  used  by  the  Draper  Laboratory 
to  determine  the  values  for  the  cross-sectional  area  A, 
the  moment  of  inertia  I,  and  the  polar  moment  of  inertia  J 
(Ref  9)  ; 


A  =• 

2irrt 

-  2ffr^ 

■  27rr^  “  Trr^ 

(4) 

m 

Tiu:  EF 

nr^ 

(5) 

XX 

yy 

UU 

J  - 

2irr^ 

-  irr^  -  21 

(6) 

rTF 

W 

The  non-structural  masses  (mirrors,  equipment  package. 
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and  solar  panels)  are  lumped  at  23  node  points;  these 
masses  and  their  locations  are  listed  as  the  C0NM2  entries 
in  the  NASTRAN  input  deck . 

The  eigenvalues  and  natural  frequencies  from  the  30 
mode  eigenvalue  analysis  performed  by  NASTRAN  are  listed  in 
Table  I.  From  these  30  modes,  the  12  modes  marked  with  an 
asterisk  were  selected  for  further  study.  These^  12  modes 
were  selected  because  they  are  the  most  observable/control¬ 
lable  or  have  the  largest  input  to  the  line  of  sight  error. 
(Ref  1:4-25)  The  B  matrix  for  these  12  modes  is  listed  in 
Appendix  H.  Since  the  sensors  and  actuators  are  collocated, 
the  C  matrix  is  the  transpose  of  the  B  matrix. 

For  this  research,  the  control  system  consists  of  21 
collocated  sensors  and  actuators  located  at  nodes  9,  10, 

11,  and  12  on  the  lower  mirror  support  truss  and  nodes  27- 
30  and  32-35  on  the  upper  mirror  support  truss.  This  con¬ 
trol  system  was  suggested  by  CSDL  in  Reference  10.  A  list¬ 
ing  of  the  node  point  and  orientation  associated  with  each 
sensor/actuator  pair  is  listed  in  Appendix  G.  Figure  3 
shows  the  orientation  and  location  of  these  pairs  on  the 
CSDL  2  model.  This  control  set  will  be  used  throughout 
this  research. 

Perturbations  to  the  Model 

Perturbations  to  the  CSDL  2  model  are  used  to  deter¬ 
mine  the  sensitivity  of  the  control  system  to  variations  in 
frequency  and  mode  shape.  These  perturbations  will  also 
aid  in  the  analysis  of  the  effects  of  system  modeling 
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TABLE  I 


Eigenvalues  and  Frequencies  for  the  First  30  Modes 


Mode  Number 

Eigenvalue (rad/ sec) 

^Frequency  (hz) 

1 

0.0 

0.0 

2 

0.0 

0.0 

3 

0.0 

0.0 

4* 

0.0 

0.0 

5* 

0.0 

0.0 

- ^ - 

0.0 

7* 

.51276 

.114 

8 

.85351 

.147 

9 

.88129 

.150 

10 

1.2120 

.175 

- IT - 

8.1889 

12* 

12.2635 

.557 

13* 

14.0320 

.596 

14 

14.9206 

.615 

15 

15.9854 

.636 

- rs - 

17* 

26.2311 

.815 

18 

26.2978 

.816 

19 

26.7716 

.823 

20 

33.0992 

.916 

Tsy^ 

7972 

22* 

53.0072 

1.159 

23 

94.9786 

1.551 

24* 

124.1206 

1.773 

25 

199.8559 

2.250 

- - 

iOO.5970 

2.254 

27 

465.4279 

3.434 

28* 

470.5175 

3.452 

29* 

627.4820 

3.987 

30 

648.1487 

4.052 
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and  manufacturing  deviations  on  the  control  of  the  structure. 
There  are  an  infinite  number  of  possible  variations  to  the 
CSDL  2  model;  this  research  concentrated  on  mass  changes 
in  the  three  mirrors  and  stiffness  changes  for  selected 
structural  elements  in  the  optical  support  truss,  The 
perturbations  were  decided  upon  using  a  combination  of 
engineering  insight  and  intuition. 

The  mass  perturbations  were  accomplished  by  adding  or 
subtracting  mass  to  various  combinations  of  the  three  mir¬ 
rors,  Because  the  VCOSS  version  of  the  CSDL  2  model  was 
designed  "to  meet  the  minimtim  requirements  for  local  sta¬ 
bility  and  frequency"  and  the  structural  elements  were 
sized  "to  prevent  buckling  or  excessive  dynamic  interaction" 
(Ref  3:17),  the  mass  changes  were  limited  to  10  percent  or 
less  of  the  mirror  mass.  This  limitation  reduces  the  im¬ 
pact  of  increased  or  decreased  mirror  mass  on  a  model  de-^ 
signed  for  specific  mirror  masses.  Larger  mass  changes 
would  result  in  larger  variations  of  frequency  and  mode 
shape,  however,  the  changed  model  would  no  longer  be  rep¬ 
resentative  of  a  real-rworld  system  and  a  major  constraint 
of  the  VCOSS  program  (realistic  design)  would  not  be  satis¬ 
fied, 

The  combinations  of  mass  changes  made  to  the  model  are 
listed  in  Table  II.  1,  5,  and  IQ  percent  changes  were  ap¬ 
plied  to  all  three  mirrors,  the  secondary  and  tertiary  mir- 


rors,  then  the  secondary  and  tertiary  mirrors  by  themselves. 
These  percentage  changes  were  added  to  and  subtracted  from 
the  original  masses.  The  primary  mirror  mass  was  not  changed 
independently  of  changes  to  the  secondary  and  tertiary  mir¬ 
ror  masses  because  the  CSDL  2  model  is  relatively  symmetric 
and  it  was  felt  that  the  results  of  changing  the  mass  of 
only  the  tertiary  mirror  would  be  similar  to  the  results  of 
changing  the  mass  of  only  the  primary  mirror ^  (Both  mirrors 
are  in  the  upper  mirror  support  truss,)  A  similar  assump¬ 
tion  was  made  with  regards  to  changing  the  masses  of  only 
the  primary  and  secondary  mirrors. 

Stiffness  changes  to  structural  elements  in  the  CSDL  2 
model  are  accomplished  by  changing  the  cross-sectional  area 
and  moments  of  inertia  of  the  epoxy-graphite  tubes,  A  per¬ 
centage  change  is  made  to  the  cross-sectional  area  and  then 
the  radius  of  the  tube  is  determined.  The  moments  of  inertia 
are  then  calculated  using  the  new  radius,  Throughout  this 
set  of  perturbations,  r/t  is  maintained  at  120.  Because 
this  model  already  uses  the  minimum  acceptable  wall  thick¬ 
nesses,  the  percentage  changes  are  only  added  to  the  values 
used  in  the  VCOSS  model  (i,e,,  the  bars  are  stiffened).  For 
example,  one  of  the  perturbations  increased  the  area  of  CBARs 
99-102  by  10  percent ,  The  new^  area  for  these  CBARs  is  (1 , 1) 
CQ,470559E-a3)  -  0, 517615E-Q3,  Using  equation  (.4),  A  «  irr^/ 
60  or  r  ■  CbOA/n)^  "  0,0994  meters.  Using  equations  C5)  and 
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(6).  I  -  0.255572E-05  and  J  -  0.511144E-05. 

A  list  of  the  values  input  to  NASTRAN  and  the  CBARs 
they  represent  is  contained  in  Table  III.  The  particular 
elements  were  selected  because  they  are  elements  that  di> 
rectly  pass  the  energy  of  the  vibrating  mirrors  to  the  rest 
of  the  optical  support  structure.  Figures  4  through  12 
show  the  position  of  the  stiffened  elements  in  the  structure. 

CBARs  44  and  45  (Figure  13)  were  stiffened  to  the  point 
of  making  them  solid  to  simulate  an  added  mass  in  the  mid¬ 
dle  of  the  optical  support  structure.  As  originally  de¬ 
signed,  the  cross-sectional  area  of  these  bars  is 
Q.104152E-Q3  m^,  r  is  0.0446  m,  and  their  length  (from  their 
node  point  coordinates)  is  8  meters.  The  material  density 
for  all  of  the  structural  elements  in  the  model  is  1720  kg/m^, 
making  the  mass  of  these  bars  1.433  kg.  This  mass  was  in¬ 
creased  by  5,  25,  and  50  kg  by  keeping  r  constant  and  increas¬ 
ing  the  wall  thickness  t.  Finally,  the  2  CBARs  were  made 
solid,  resulting  in  a  mass  of  85.988  kg  per  bar.  Table  IV 
lists  the  NASTRAN  input  data  for  these  stiffness/mass  changes. 
The  values  for  the  25  kg  increase,  for  example,  were  derived 
as  follows: 

mass  ■  (area) Clength) Cdensity) 

Area  A  -  mass  -  1.433  +  25  -  .001912  m^ 


mass 

Length  x  density 


A  ■  2Trrt 


I 


TABLE  III 

Stiffness  Perturbations  to  the  CSDL  2  Model 


Set 


1 

2 

3 

4 

5 

6 

7 

8 

9: 


Percent  changes  in  Area  added  to  CBARs 


99-102 

99-102  &  87-97  (odd) 
90  and  94 
90 

123-126 

123-124 

76,78,80,81,83,85 
76,78,80,81,83,85  & 

87- 97  (odd) 

88- 98  (even) 


Perturbation^% 
Set  Number  y 

Increase 
n  A 

Old  A 

New  A 

New  ljQj.“lyy 

New  J 

1 

10 

.470559E-03 

.517615E-03 

.255572E-05 

.511144E-05 

50 

II 

.705838E-03 

.475573E-05 

.951507E-05 

2 

10 

It 

.517615E-03 

.255572E-05 

.511144E-05 

.259503E-03 

.285453E-03 

.778109E-06 

.155622E-05 

50 

.470559E-03 

.705838E-03 

.475573E-05 

.951507E-05 

.259503E-03 

.389255E-03 

.144690E-05 

.289381E-05 

3 

10 

.566323E-03 

.622955E-03 

.370582E-05 

.741165E-05 

50 

ft 

.849485E-03 

.689101E-05 

.137820E-04 

4 

10 

II 

.622955E-03 

.370582E-05 

.741165E-05 

50 

II 

.849485E-03 

.689101E-05 

.137820E-04 

5 

10 

.716806E-04 

.788487E-04 

.593691E-07 

.118738E-06 

50 

II 

.107521E-03 

.110397E-06 

.220794E-06 

6 

10 

II 

.788487E-04 

.593691E-07 

.118738E-06 

50 

It 

.107521E-03 

.110397E-06 

.220794E-06 

7 

10 

.259503E-03 

.285453E-03 

.778109E-06 

.155622E-05 

50 

II 

.389255E-03 

.144690E-05 

.289381E-05 

8 

10 

II 

.285453E-03 

.778109E-06 

.155622E-05 

50 

ft 

.389255E-03 

.144690E-05 

.289381E-05 

9  (88&96) 

10 

.325124E-03 

.357636E-03 

.122139E-05 

.244278E-05 

(90&94) 

.566323E-03 

.622955E-03 

.370582E-05 

:741146E-05 

(92&98) 

.348640E-03 

.383504E-03 

.140447E-05 

.280893E-05 

(88&96) 

50 

.325124E-03 

.487686E-03 

.227118E-05 

.454236E-05 

(90&94) 

.566323E-03 

.849485E-03 

.689101E-05 

.137820E-04 

(92&98) 

.348640E-03 

.522960E-03 

.261161E-05 

.522322E-04 
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Figure  4.  Perturbation  Set  1 
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Figure  5.  Perturbation  Set  2 
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I 


Figure  8.  Perturbation  Set  5 


i 

! 


Figure  9,  Perturbation  Set  6 


Figure  12.  Perturbation  Set  9 


CBARs  44  and  45 
(Midsection  Mass  Increases) 


t  =  A  -  ,001912  =  .006855  tn 

2Trr  Ztt  (  0446) 

r/t  =  .0446  =  6.5061 

700^55 

I  =  Trr^  =  Tr(.0446)^  =  .191062E-05 

FTE  ~Tr5mr 

and 

J  -  21  =  .382125E-05 

•  s 

The  analysis  of  changes  in  the  damping  factor  will  use 
data  from  the  three  perturbation  groups  just  discussed,  The 
appropriate  changes  are  made  to  the  NASTRAN  input  deck  and 
the  output  is  processed  through  DOFB34.  The  equation  used 
to  calculate  the  damping  factor,  as  well  as  those  used  to 
find  the  B  and  C  matrices,  are  the  subject  of  the  next  chap¬ 
ter  of  this  thesis. 
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III.  Theory 


B  and  C  Matrices 

NASTRAN  is  rvin  to  determine  the  modal  frequencies  and 
mode  shape  values  at  node  points  on  the  structure.  The  re¬ 
sults  of  NASTRAN  are  the  matrix,  which  is  a  row  matrix 
containing  the  eigenvectors  (mode  shapes)  of  the  structure, 
and  a  listing  of  the  modal  frequencies. 

The  equations  of  motion  for  a  vibrating  structure  such 
as  CSDL  2  are 

Mq  +  Eq  +  Kq  =  Du  (7) 

Where  q  is  an  n-vector  of  generalized  coordinates  that  com¬ 
pletely  describe  the  motion  of  the  system,  M  an  n  x  n  sym¬ 
metric  mass  matrix,  K  an  n  x  n  symmetric  stiffness  matrix, 

E  an  n  X  n  symmetric  damping  matrix,  u  an  m-vector  of  in¬ 
puts  (disturbances) ,  and  D  an  n  x  m  matrix  of  actuator  co¬ 
efficients.  The  vector  Du  therefore  represents  the  input 
force  vector.  (Ref  6:1)  Equation  (7)  can  be  rewritten  in 
model  coordinates,  3,  where  q  *  ij)n.  The  mode  shape  (repre¬ 
sented  by  <j>)  is  unique,  but  the  amplitude  is  not.  By  re¬ 
writing  equation  (.7)  in  modal  coordinates,  the  elements  of 
the  natural  modes  are  made  unique  (i.e.,  they  are  normalized). 
In  this  case,  this  was  accomplished  by  setting  (f)TH(|>  equal  to 
I.  (Ref  5:141)  Equation  (7)  is  now  written 


Where  u  is  the  natural  (undamped)  modal  frequency  and  t;  is 
the  damping  factor.  The  matrix  |^2c(ijJ  is  therefore  an  n  x  n 
diagonal  damping  matrix  and  is  an  n  x  n  diagonal  matrix 

of  the  eigenvalue  of  equation  C7) .  Because  -  I,  the 

modal  matrix  4)  is  such  that  and  -  [zcuj  . 

(Ref  5:141)  Equation  (8)  can  be  written  in  the  state  vector 
form 

X  “  Ax  +  Bu  (.9) 


Where 


nT. 


and 


B 


The  general  output  equation  is  of  the  form 

Z  “  Cp£  +  Cvi 


(10) 


assuming  that  both  position  and  velocity  sensors  are  being 
used.  In  terms  of  x, 


Where 


C  »  {Cp4),  C^ij) } 


Cp  is  a  matrix  corresponding  to  the  direction  cosines  of  the 
position  sensors  and  corresponds  to  the  direction  cosines 
of  the  velocity  C^ate)  sensors.  The  designated  control  set 
for  this  research  uses  only  collocated  rate  sensors  and  force 
actuators  (no  position  sensors)  .  Thus ,  Cp  -  0  and  -  C 
where  the  rows  of  the  B  matrix  physically  correspond  to  the 
amplitude  of  each  mode  at  the  various  actuator  locations  and 
the  columns  of  C  represent  the  amplitude  of  each  mode  at  the 
various  sensor  locations.  Negating  the  zero  half  of  the  B 
matrix  (from  equation  (9))  gives 

B  •  (12) 

Therefore,  an  element  of  B  can  be  written 

■  V,!  '““j  Vl.t  ““j  V2.1 

where  the  cosine  terms  are  the  direction  cosines  of  the  D 
matrix.  Equation  (13)  is  used  in  the  ZETA  program  to  calcu¬ 
late  the  B  and  matrices. 

Damping  Factor 

Program  D0FB34  calculates  two  eigenvector  (complex 


conjugate)  pairs  per  mode,  each  having  a  real  value  and  an 
imaginary  value.  The  magnitude  of  the  real  value  is  the 


product  of  the  damping  factor  and  the  natural  frequency; 
the  magnitude  of  the  imaginary  value  is  the  damped  frequency 
‘**d.  In  Figure  14,  one  of  the  eigenvalues  is  plotted  in  the 
s-plane.  From  this  depiction 

tan  n  ■ 

or  n.  "  tan~^  u*  ■  cos”^? 

and  C  »  cos  tan~^  oj'  (14) 

P" 

Equation  (14)  is  used  to  calculate  the  damping  factor 
for  each  mode  that  is  analyzed  in  this  research.  One  of  the 
inputs  to  D0FB34  is  the  damping  factor  of  the  uncontrolled 
structiire  that  occurs  as  a  result  of  the  structure's  resis¬ 
tance  to  bending  and  twisting.  For  this  research,  a  value 
of  0.01  was  assigned.  The  purpose  of  a  control  system  is 
to  actively  danqjen  vibrations  in  the  structure.  Thus,  an 
increase  in  the  damping  factor  Ci>e.,  a  value  greater  than 
0.01)  means  that  the  control  system  is  actively  damping  that 
particular  mode  to  some  extent.  If  the  value  stays  at  or 
very  near  .01,  the  control  set  is  having  little  or  no  effect 
on  that  particular  mode. 

Using  eqxiatlons  developed  in  this  chapter,  data  points 
for  determining  the  effects  of  perturbations  on  the  damping 
factor  can  be  generated.  The  next  chapter  presents  this 


complex  plane  (Ref  11:116) 


! 

i 


data  and  discusses  the  impact  that  the  control  system  ori¬ 
entation  and  the  modeling  of  the  structure  have  on  the  find¬ 
ings. 
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IV.  Results 


In  this  chapter,  the  mode  shapes  for  the  twelve  modes 
are  discussed  and  the  effect  of  the  perturbations  on  the 
damping  factor  for  each  mode  are  presented. 

Mode  Shapes 

Appendix  I  (Figures  15  to  24)  contains  views  of  the 
twelve  mode  shapes  from  the  vmperturbed  model  under  consid¬ 
eration  in  this  research.  These  plots  were  made  using 
GCSNAST,  a  program  that  uses  the  NASTRAN  bulk  data  deck 
and  a  deflection  file  that  is  derived  from  the  NASTRAN- 
produced  modal  frequencies  and  mode  shapes.  (See  Appendix 
K  for  GCSNAST  procedures.)  In  these  plots,  the  dashed 
lines  represent  the  deformed  structure  and  the  solid  lines 
represent  the  original  structure.  The  deformations  have 
been  magnified  by  a  factor  of  ten  for  clarity.  A  brief 
description  of  the  twelve  modes  follows. 

Modes  4,  5,  and  6  (Figure  15)  are  rigid  body  rotation 
modes  encompasing  rigid  body  rotation  around  each  of  the 
three  axes. 

Mode  7  (Figure  16)  is  a  mode  where  the  solar  panels 
are  deforming  in  the  xy  and  xz  planes  and  the  structure 
rotates  in  the  xz  plane  about  a  point  near  the  middle  of 
the  equipment  section.  There  is  no  noticeable  deformation 
in  the  structure. 
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Mode  12  (Figure  17)  involves  rotation  and  bending  in 
the  xy  and  xz  planes  and  a  considerable  deformation  of 
the  solar  panels .  The  point  of  rotation  is  in  the  upper 
part  of  the  optical  support  structure. 

Mode  13  CFigxire  18)  has  rotation  in  the  yz  plane  and 
deformation  of  the  solar  panels  in  the  xy  and  xz  planes. 

The  rotation  is  about  a  point  in  the  upper  mirror  support 
truss . 

Mode  17  (Figure  15)  involves  a  slight  rotation  in  the 
yz  plane  and  deformation  of  the  solar  panels  in  the  xy  and 
xz  planes.  The  structure  rotates  about  a  point  near  the 
outside  edge  of  the  primary  mirror. 

In  Mode  21  (Figure  20)  there  is  twisting  and  solar 
panel  deformation  in  the  xy  plane  and  deformation  of  all 
three  components  of  the  optical  support  structure  in  all 
three  planes. 

Mode  22  (Figure  21)  involves  twisting  and  solar  panel 
deformation  in  the  xz  plane. 

Mode  24  (Figure  22)  has  the  most  radical  deformation 
of  the  twelve  modes.  There  is  twisting  in  the  xy  and  yz 
planes  and  deformation  of  the  optical  support  structure 
and  solar  panels  in  all  three  planes. 

Except  for  the  solar  panels,  the  deflections  in  mode 
28  (Figure  23)  are  quite  small.  There  is  some  deformation 
in  the  upper  and  lower  mirror  support  trusses  in  the  xy  and 
yz  planes. 

Mode  29  (Figiire  24)  also  has  very  small  deformations 
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in  the  optical  support  structure.  The  metering  truss  is 
deformed  slightly  in  the  yz  plane  and  the  solar  panels  are 
deformed  in  the  xy  and  xz  planes. 

In  the  CSDL  2  model,  the  equipment  section  is  connected 
to  the  optical  support  structxire  by  three  springs  at  nodes 
3,  4,  and  6.  These  springs  totally  isolate  the  equipment 
section  from  the  optical  support  structure  and,  therefore, 
the  energy  from  the  vibrating  solar  panels  Ccantilevered 
from  the  equipment  section),  is  not  passed  upward  to  the 
optical  support  structure. 

Damping  Factor 

Appendix  J  lists  the  damping  factors  calculated  during 
this  research.  The  damping  factors  associated  with  modes 
12,  13,  17,  and  24  of  the  perturbed  models  are  portrayed 
in  Figures  25  to  36.  In  these  figures,  the  vertical  scale 
is  the  damping  factor  in  percent;  the  horizontal  scale  is 
the  percent  change  in  mass  (for  the  mirror  mass  pertur¬ 
bations)  ,  the  percent  increase  in  cross-sectional  area 
(for  the  stiffness  perturbations) ,  or  the  amount  of  mass 
that  was  added  to  the  middle  of  the  optical  support  struc¬ 
ture.  The  horizontal  scale  is  not  linear  and  the  lines 
connecting  the  various  data  points  are  there  only  to  link 
the  data  points  for  a  particular  mode.  Data  points  on  the 
vertical  axis  are  those  for  the  unperturbed  structure. 

The  other  eight  of  the  original  twelve  modes  are  not  plotted 
because  their  damping  factors  did  not  change  or  did  not 
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contribute  to  the  analysis. 


Modes  4,  5,  and  6  are  rigid  body  modes  and  have  a 
damping  factor  of  zero.  The  eigenvalue  pairs  for  these 
three  modes  lie  on  the  imaginary  axis  (See  Figure  14)  and 
the  response  has  sustained  oscillations  (Ref  11:242).  The 
designated  control  set  for  this  analysis  uses  only  rate 
(velocity)  sensors.  Since  the  natural  frequency  for  a 
rigid  body  mode  is  zero  (the  structure  only  translates  or 
rotates,  it  does  not  deform)  the  w  terms  in  the  A  matrix 
(eqxiation  (9))  are  zero  and  A  is  therefore  zero.  The  sen¬ 
sors  are  sensing  the  motion  but  there  is  no  control  system 
input  due  to  A  being  zero.  The  rigid  body  modes  will  not 
be  controllable  unless  position  sensors  are  added  to  the 
control  set. 

The  eigenvalue  pairs  for  mode  7  lie  on  the  real  axis 
but  are  not  identical,  hence  the  damping  factor  is  greater 
than  one,  the  mode  is  overdamped,  and  the  oscillation  is 
damped  out  in  less  than  one  cycle.  The  control  system 
(Figure  3)  is  oriented  such  that  there  are  three  sensors/ 
act\iator  pairs  that  sense  motion  in  the  x-direction  and 
twelve  in  the  z-direction,  all  of  which  are  affected  by 
this  mode  shape.  This  causes  the  overdamping  in  this  mode. 
The  overdamping  remains  throughout  the  research,  so  mode  7 
will  not  be  analyzed  further. 

Throughout  the  perturbations,  the  damping  factors  for 
modes  22,  28,  and  29  have  remained  at  .01,  which  is  the  value 
assumed  to  exist  in  the  basic  uncontrolled  structure.  Modes 


28  and  29  show  very  little  deformation  in  the  optical  sup¬ 
port  structure.  Because  of  the  small  deformations,  the 
sensors  do  not  deflect  very  much.  This  limits  their  input 
to  the  control  system  processors  and,  hence,  limits  the 
control  systems  effectiveness  in  damping  these  modes.  Mode 
22  is  primarily  deformation  In  the  z-direction,  where  there 
are  twelve  sensor/actuator  pairs  oriented  to  sense  motion. 
But,  because  of  the  small  deformations  and  the  location  of 
the  sensors ,  the  control  system  can  not  effectively  damp 
the  vibration.  These  three  modes  C22,  28,  and  29)  lie  in 
iincontro liable  modal  space  and  are  not  controllable  given 
the  orientation  of  this  particular  control  system.  That  is, 
the  location  and/or  orientation  of  the  sensors  has  to  be 
changed  such  that  they  can  more  completely  sense  the  motion 
of  the  vibration.  Except  for  this  observation,  these  modes 
present  no  other  useful  data  for  this  research  and  will  not 
be  analyzed  further. 

The  damping  factor  for  mode  21  varies  by  only  .00007 
in  this  research  (.01015  to  .01022)  and  will  be  placed  in 
the  same  category  as  modes  22,  28,  ard  29.  Thus,  the 
analysis  will  concentrate  on  modes  12,  13,  17,  and  24  and 
focus  on  the  changes  in  the  damping  factor  with  respect  to 
the  physical  changes  (perturbations)  in  the  CSDL  2  model. 
Changes  in  the  damping  factor  of  less  than  .001  (10%  of 
the  passive  damping  factor,  .011  will  be  assumed  to  be 
insignificant. 


In  some  cases,  there  is  a  definite  trend  in  the  damp¬ 
ing  factor  data  as  the  magnitude  of  the  perturbation  is 
increased  or  decreased  CFigure  25,  Mode  13;  Figure  31, 
Mode  13;  Figiare  32;  Figure  33,  Mode  13;  Figure  34,)  Con¬ 
versely,  there  are  many  cases  where  the  change  in  the 
damping  factor  does  not  follow  a  definite  pattern  CFigure 
25,  Mode  12;  Figure  26,  Mode  17;  Figure  27  to  30;  Figure 
31,  Modes  12  and  17;  Figures  35  and  36).  After  looking 
at  how  the  damping  of  each  mode  is  affected  by  the  three 
basic  types  of  perturbations,  there  will  be  a  discussion 
on  the  modeling  of  the  structure,  the  orientation  of  the 
control  system,  and  their  impact  on  the  results. 

Mode  12 :  Due  to  the  sensor/actuator  configuration, 
there  is  only  one  sensor  (at  node  11)  that  senses  the 
motion  of  the  lower  part  of  the  structure  in  the  x- 
direction.  There  is  some  twisting  in  the  xy  plane.  This 
twisting  is  most  notable  in  the  upper  mirror  support  truss, 
and,  from  Figure  17,  it  is  noted  that  the  tertiary  mirror 
is  deflected  more  than  the  primary  mirror.  There  are 
sensors  at  node  34  (near  the  tertiary  mirror)  that  sense 
X  and  y  motion,  and  node  34  is  at  the  point  of  maximum 
deflection.  Both  deflections  that  occur  in  mode  12  are, 
therefore,  sensed  by  the  control  set  and  it  is  anticipated 
that  there  will  be  some  degree  of  active  damping  in  this 
mode. 

The  actual  effects  of  mirror  mass  changes  on  this  mode 
are  portrayed  in  Figtires  25  and  27.  The  only  combination 
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whose  damping  factor  was  unchanged  or  better  in  every  case 
is  the  addition  of  mass  to  the  secondary  mirror.  It  should 
be  noted,  however,  that  the  +57.  value  is  greater  than  the 
+10%  value.  Addition  of  mass  to  all  three  mirrors  and  to 
both,  the  secondary  and  tertiary  mirrors  produced  lower 
damping  factors  in  every  case  except  for  the  +10%  value 
for  the  two  mirrors.  The  other  perturbations  produced  no 
trends  in  that  the  damping  factors  were  not  consistently 
higher  or  lower  than  the  vinperturbed  value  as  the  masses 
were  increased.  The  changes  to  the  secondary  and  tertiary 
mirrors  by  themselves  produced  results  that  were  almost 
identical,  hut  they  do  not  show  any  consistent  trend. 

The  results  for  the  stiffness  changes  are  shown  in 
Figures  29.  to  34.  Perturbation  sets  2,  5,  6,  and  9  exhibit 
consistent  trends  (The  damping  factor  stays  the  same  or 
gets  better/worse  for  both  data  points.).  Set  2  stiffens 
the  upper  six  vertical  members  in  the  metering  truss 
(Figure  5) .  Sets  5  and  6  (figures  8  and  9)  stiffen  all 
four,  then  two,  of  the  upper  diagonal  members  that  link 
the  outer  edges  of  the  primary  and  tertiary  mirrors  to  the 
top  of  the  structure.  It  is  noted  that  set  6,  which 
stiffens  just  the  two  members  supporting  the  tertiary 
mirror,  produces  the  same  results  as  set  5.  This  is 
consistent  with  the  observation  that  the  tertiary  mirror 
is  most  affected  by  the  twisting  motion.  The  data  from 
these  two  sets  implies  that  it  is  the  stiffening  of  the 
tertiary  mirror  supports  that  causes  the  changes  in  the 
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damping  factor.  Set  9.  (.Fig^lre  12)  stiffens  the  six  diag¬ 
onal  supports  in  the  upper  portion  of  the  metering  truss. 

This  could  increase  the  amoxint  of  the  twisting  motion  that 
is  passed  downward.  There  are  no  sensor/actuator  pairs  to 
damp  this  motion  in  the  middle  of  the  structTire,  so  the 
dancing  factor  would  decrease.  Sets  1,  3,  4,  7,  and  8  do 
not  give  consistent  results. 

The  results  from  the  midsection  mass  changes  are 
shown  in  Figure  35.  The  process  of  adding  mass  has  stiff¬ 
ened  the  two  elements  shown  in  Figure  13,  so  it  is  possible 
that  this  could  increase  the  resistance  to  twisting  in  the 
z-direction.  In  this  mode,  the  damping  increases  after 
approximately  5Q  kg  have  been  added  to  each  element,  in¬ 
dicating  that  at  this  point,  the  mot ion /inertia  in  the 
middle  of  the  structure  is  helping  to  damp  the  rotation. 

Mode  13 ;  The  major  component  in  this  mode  is  rota¬ 
tion  in  the  yz  plane.  This  motion  is  large  enough  to 
result  in  significant  displacement  in  the  y  and  z  directions. 
The  primary  mirror  is  deflected  more  than  the  tertiary 
mirror.  The  structure  is  also  elongated  in  the  z-direction. 
There  are  twice  as  many  sensor/actuator  pairs  working  in 
the  y-direction  as  there  are  in  the  x-dlrection  and  twelve 
pairs  oriented  in  the  z-direction.  Thus,  it  is  expected 
that  the  damping  factors  for  this  mode  will  be  larger  than 
for  mode  12. 

The  results  of  the  mass  changes  for  mode  13  are  por¬ 
trayed  in  Figures  25  and  27.  The  data  for  all  three  mirrors 


depicts  a  consistent  trend  xintil  the  +10%  value.  When 
comparing  the  data  from  the  two  perturbations  portrayed  in 
Figure  25 ,  it  is  noted  that  balanced  increases  in  the  upper 
mirror  support  truss  adding  eqiial  percentages  of  mass 

to  both  the  primary  and  tertiary  mirrors)  results  in  better 
damping  than  mbalanced  increases.  The  data  for  increases 
in  the  tertiary  mirror  by  itself  CFigure  27)  concurs  with 
this  finding.  It  is  also  noted  that  the  damping  factors  for 
mode  13  are  considerably  higher  than  for  mode  12,  as  expected. 

Increasing  stiffness  resulted  in  lower  damping  for  every 
stiffness  perturbation  except  for  the  +10%  change  in  set  8, 
where  there  was  a  barely  significant  increase  of  .Q014SI. 

All  of  the  sets  except  1  and  8  showed  consistent  results. 
Increasing  the  stiffness  increases  the  frequency  of  the 
vibration,  and  this  research  has  shown  that  lower  frequen¬ 
cies,  with  their  larger  deformations,  are  generally  more 
easily  controlled.  The  results  from  the  mode  13  stiffness 
perturbations  confirm  this  in  almost  all  cases. 

The  results  from  the  midsection  mass  increases 
(Figure  35)  are  consistent  except  for  the  +50  kg  data 
point.  Increasing  mass  in  the  middle  of  the  structure 
generally  lowers  the  damping. 

Mode  17 ;  In  this  mode,  there  is  a  small  rotation  in 
the  yz  plane  about  a  point  near  the  primary  mirror.  The 
deflections  in  the  mirror  support  trusses  are  quite  small 
and  there  are  only  two  sensor/actuator  pairs  in  the  lower 
part  of  the  structure  that  can  sense  this  motion.  The 


displacements  in  the  z-direction  are  very  small,  so  the 
twelve  pairs  that  sense  motion  in  the  z-direction  will  not 
be  displaced  very  much  and  are  not  expected  to  have  much 
effect  on  the  damping  of  this  mode.  Damping  is  expected 
to  be  lower  than  for  modes  12  and  13 . 

Data  for  the  mirror  mass  changes  is  shown  in  Figures 
26  and  28.  In  all  cases,  the  damping  factor  increased  or 
remained  the  same.  There  is  not  a  consistent  trend  as  the 
mirror  masses  were  varied  from  -10%  to  +10%  of  their  original 
values.  Every  perturbation  group  (i.e.,  all  mirrors,  ter¬ 
tiary  only,  etc.)  has  at  least  two  data  points  that  are  in¬ 
consistent  with  the  others. 

The  stiffness  changes  (Figure  29  to  34)  produced  more 
consistent  results  and  point  out  two  perturbations  that  had 
a  very  detrimental  effect  on  the  damping  of  this  mode.  The 
damping  factors  were  xonchanged  or  higher  except  for  those 
two  cases.  Sets  2,  5,  6,  and  9  showed  increasing  damping 
as  stiffnesses  were  increased.  Sets  1,  3,  and  4  had  more 
damping  at  the  10%  increase  in  area  than  at  the  507.  point . 

The  two  detrimental  cases,  the  107.  increases  for  sets  7  and 
8,  showed  almost  no  active  d2uiiping  at  all.  Both  of  these 
sets  increased  the  area  of  the  lower  six  vertical  members 
of  the  metering  truss;  set  8  also  increased  the  area  of  the 
corresponding  vertical  members  in  the  upper  part  of  the 
metering  truss.  It  would  appear  then  that  the  structure 
is  very  sensitive  to  stiffness  changes  in  the  lower  ver¬ 
tical  members.  When  the  areas  of  the  members  in  these 
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two  perturbation  sets  were  increased  by  507o,  active  damp¬ 
ing  returned  and  the  damping  factor  was  greater  than  for 
the  tinperturbed  structure. 

One  data  point  in  the  midsection  mass  increases 
(Figure  36)  showed  a  drastic  decrease  in  damping.  Except 
for  this  point  (25  kg  increase) ,  the  damping  was  the  same 
or  better  than  that  of  the  mperturbed  case. 

Mode  24:  This  mode  is  ^characterized  by  large  deforma¬ 
tions  in  all  three  planes.  Thus,  all  21  sensor/actuator 
pairs  get  involved  and  there  is  some  degree  of  active 
damping  even  though  the  frequency  (1.773  hz)  is  higher  than 
the  other  modes  where  active  damping  was  subsequently 
employed.  The  damping  factors  for  all  three  perturbation 
groups  are  lumped  very  closely  around  a  value  of  .012, 
with  no  data  points  varying  by  more  than  the  .001  value 
that  has  been  defined  as  the  significant  amotint  for  this 
research.  It  can  be  stated,  then,  that  mode  24  is  actively 
damped  and  that  the  perturbations  had  no  significant  effect 
on  the  damping  factor.  It  is  possible  that  more  sensor/ 
actuator  pairs  oriented  in  the  x  and  y  directions  would 
increase  the  damping  for  this  mode. 

This  chapter  has  discussed  the  mode  shapes  and  the 
changes  in  the  damping  factor  as  the  structure  was  perturbed. 
The  next  chapter  will  discuss  the  possible  causes  for  the 
lack  of  consistent  results,  summarize  the  findings,  and 
present  recommendations  for  further  research  to  remove  or 
lessen  the  uncertainty  in  these  results. 
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CHANGE  IN  MIRROR  MASS ,  7. 

Figure  25.  Mirror  Mass  Changes 
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INCREASE  IN  CROSS-SECTIONAL  AREA,  7. 


Figure  29.  Stiffness  Changes,  Set  1 
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Figure  31.  Stiffness  Changes,  Sets  3  and  4 


Figure  32.  Stiffness  Changes,  Sets  5  and  6 
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Figure  33.  Stiffness  Changes,  Sets  7  and  8 
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Figure  34.  Stiffness  Changes,  Set  9 


MASS  ADDED,  kg 

Flgtire  35,  Mid-section  Mass  Increases 


K^SS  ADDED,  kg 

Figure  36.  Mid-section  Mass  Increases 


V.  Conclusions  and  Recommendations 


Discussion  on  Inconsistencies  in  the  Data 

Due  to  the  inconsistency  of  portions  of  the  data,  the 
modeling  of  the  CSDL  2  structure  and  the  particular  VCOSS 
derivation  were  studied.  In  Revision  1  of  the  CSDL  2 
model,  the  cross-members  that  go  from  node  points  28  to 
30  and  from  32  to  33  are  strengthened  to  support  the 
masses  of  the  primary  and  tertiary  mirrors,  respectively. 
(This  research  used  the  revision  1  model).  The  cross- 
section  of  the  stiffened  mirror  support  beams  is  depicted 
in  Figure  37.  They  are  tubular  trusses  and  were  designed 
to  prevent  line-of-sight  errors  due  to  bending  at  frequen¬ 
cies  less  than  40  hz.  (Ref  3:9)  The  cross  members  and 
exterior  members  of  these  trusses  do  not  run  the  full 
length  of  the  truss  and  are  assumed  to  have  no  effect  on 
its  bending  stiffness.  The  mass  per  mit  length  is  set 
at  1.5  times  that  of  the  four  tubes  by  themselves  to 
account  for  the  cross  members.  In  the  NASTRAN  input  deck, 
this  tubular  truss  is  modeled  as  a  circular  beam  with  the 
same  cross-sectional  area  as  the  four  tubes  combined.  The 
moments  of  inertia  for  the  equivalent  circular  beam  are  the 
same  as  the  summed  moments  of  inertia  for  the  individual 
tubes  aroxmd  the  center  of  the  truss.  The  Inertia  of  the 
cross  members  and  exterior  members  has  not  been  included. 

The  asaxjmption  has  been  made  through  this  modeling  technique 
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Tube  Properties : 

Radius  =  0 . 11m 

t  -  0.0022m 
Area  =■  0.00152m^ 


Section  Properties: 

Area  =  4A^  ,  =■  0.00608m‘^ 

tube  /, 

I  =  0.00152m^ 

J  -  0.00304m^ 

Mass/L  =>  15.69  kg/m 


that  the  equivalent  circular  beam  will  bend  and  twist  in 
the  same  manner  as  the  tubular  truss.  This  is  a  very 
questionable  structural  modeling  technique  and  could  be 
a  factor  in  the  inconsistency  of  the  damping  factor  data 
from  this  research. 

The  approach  used  in  tha  design  of  the  VCOSS  version 
of  the  CSDL  2  model  was  to  reduce  the  size  of  the  struc- 
ttnral  elements  to  the  minimum  values  that  would  still  pre¬ 
vent  structural  failure.  The  structure  is  therefore  very 
flexible  and  has  very  low  natural  frequencies.  (Ref  3:17) 
Resizing  these  elements  downward  from  the  basic  model 
Crevision  11  lowered  the  total  mass  by  373  kg  to  8963  kg, 
a  reduction  of  slightly  over  4%.  If  an  LSS  such  as  this 
telescope  was  to  be  launched  into  space  on  the  shuttle, 
the  LSS  mass  of  8963  kg  would  be  small  compared  to  the 
shuttle  capacity  of  around  27000  kg.  A  savings  of  373  kg 
could  be  insignificant  due  to  the  fact  that  the  tubular 
elements  would  occupy  a  large  volume  in  the  payload  bay 
without  approaching  the  constraint  on  total  payload  weight. 
The  VCOSS  structure,  then,  is  perhaps  unrealistic  and  does 
not  model  the  characteristics  of  the  type  of  LSS  that  will 
be  erected  in  the  future.  This  research  has  indicated 
that  the  structure,  as  it  exists  in  the  VCOSS  configuration, 
does  not  exhibit  much  consistency  in  its  response  to  small 
changes  in  mass  and  stiffness.  This  could  indicate  serious 
problems  in  its  basic  design,  and,  since  the  goals  of  the 
VCOSS  program  include  realistic  design  and  control  hardware 


specification,  the  effectiveness  of  the  VCOSS  program  could 
be  degraded. 

The  control  system  of  21  sensor/actuator  pairs  emphasizes 
the  focusing  of  the  telescope  Crather  than  the  pointing  of  it) 
in  that  twelve  pairs  are  oriented  to  sense  motion  in  the  z- 
direction.  G^isplacements  in  the  z -direction  would  cause 
the  mirrors  to  deflect  and  result  in  an  unfocused  image,) 

There  are  three  sensor/actuator  pairs  for  the  x-direction 
and  six  pairs  for  y.  The  analysis  of  modes  21,  24,  28,  and 
29  has  shown  that  it  is  difficult  to  actively  damp  these 
modes  because  of  their  higher  frequencies  and  smaller  dis¬ 
placements.  Mode  24,  however,  is  actively  damped;  it 
involves  motion  of  many  more  sensor/actuator  pairs  than  do 
the  other  three  uncontrollable  modes  due  to  the  severity 
of  its  deformation.  It  is  possible,  then,  that  if  more 
sensor/actuator  pairs  were  oriented  in  the  x  and  y  direc¬ 
tions,  modes  21,  28,  and  29  could  be  actively  damped.  It 
is  important  to  damp  these  modes,  because,  of  the  twelve 
modes  studied  in  this  research,  modes  21,  24,  and  29  have 
the  largest  impact  on  the  line-of-sight  problem  that  is 
critical  in  a  telescope.  Dtiring  a  Lockheed  Missiles  and 
Space  Company,  Inc.,  control  design  effort  using  CSDL  2, 
a  maximum  llne-of-slght  error  of  5Q  nanoradians  Cnrad)  was 
specified.  It  was  found  that  modes  21,  24,  and  29  con¬ 
tributed  633,  815,  and  630  nrad,  respectively,  to  the  line- 
of-sight  error  during  this  analysis.  (Kef  1:4-23)  These 
values  are  significantly  higher  than  the  allowable  limit 
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and  emphasize  the  need  to  damp  these  modes. 

The  methodology  used  in  this  research  assxnned  that  the 
perturbations  to  the  CSDL  2  model  do  not  significantly 
change  the  mode  shapes  or  frequencies  from  their  values 
in  the  vmperturbed  model.  This  is  important  because 
NASTRAN  orders  its  output  as  a  function  of  the  eigenvalue 
of  the  mode,  starting  with,  the  lowest  value.  (The  eigen¬ 
value  is  the  square  of  the  natural  frequency.)  The  rest 
of  the  data  reduction  assvuaes  that  the  frequencies  do  not 
change  enough  to  cause  the  mode  shapes  to  swap  positions 
in  the  output.  For  example,  the  natural  frequency  for 
mode  7  is  .145  hz;  for  mode  8  it  is  .263  hz.  If  the  per¬ 
turbations  to  the  model  caused  the  frequency  associated 
with  mode  shape  8  to  be  less  than  the  frequency  of  mode  7, 
the  output  from  NASTRAN  would  have  the  two  modes  reversed 
and  the  analysis  of  the  results  would  be  comparing  the 
damping  factor  for  two  different  mode  shapes.  This  com¬ 
parison  would  be  invalid.  To  verify  that  the  results  of 
this  research  do  not  reflect  the  above  problem,  three 
additional  sets  of  GCSNAST  mode  shape  plots  (similar  to 
those  in  Figures  15  to  24)  were  made.  These  three  sets 
are  plots  of  the  mode  shapes  from  three  of  the  perturbations 
whose  damping  factors  varied  significantly  from  the  rest  of 
the  data  points  in  the  particular  perturbation  group.  These 
points  seemed  to  be  the  most  likely  to  be  affected  by  a 
swap  in  the  order  of  mode  shapes.  They  are  (1)  the  addition 
of  1Q%  mass  to  the  secondary  mirror,  (2)  the  set  that  made 
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CBARs  44  and  45  Cnear  the  middle  of  the  structure)  solid, 
and  (3)  the  507o  increase  in  area  for  stiffness  perturbation 
set  1  (CBARs  99-102  in  the  upper  part  of  the  metering  truss) . 
Inspection  of  the  GCSNAST  plots  showed  that,  even  when  the 
mode  shapes  are  magnified  by  a  factor  of  10,  there  are  no 
noticeable  differences  in  mode  shapes  between  these  three 
perturbations  and  the  vinperturbed  structure.  (There  must 
be  some,  otherwise  the  damping  factor  would  not  change.) 

The  important  point  to  be  made  is  that  the  overall  shapes 
did  not  change,  indicating,  for  example,  that  the  mode 
shape  associated  with  the  twelvth  eigenvalue  in  the  NASTRAN 
output  remained  essentially  the  same  throughout  the  research 
and  that  the  analysis  is  not  comparing  the  damping  for  two 
different  mode  shapes.  The  grouping  of  the  data  (i.e.,  the 
damping  factors  for  mode  12  are  all  near  .185,  mode  13  is 
around  .22,  etc.)  supports  this  conclusion.  The  assumption 
remains,  however,  that  the  mode  shapes  for  modes  that  have 
frequencies  very  close  together  are  different. 

Conclusions 

This  research  has  shown  that  the  damping  in  a  LSS  can 
be  affected  by  very  small  changes  in  mass  and  stiffness.  A 
structure  in  space  could  be  "fine  tuned"  to  increase  the 
damping  for  a  particular  mode.  It  was  noted,  however, 
that  not  all  modes  are  affected  in  the  same  way.  One  mode 
may  show  increased  damping  while  another  modes  damping  may 
decrease.  Any  attempt  to  increase  damping,  therefore,  must 


consider  all  of  the  modes  that  are  to  be  controlled. 

The  amount  of  active  damping  in  the  CSDL  2  model  is 
dependent  upon  the  mode  shape,  the  natural  frequency,  and 
the  orientation  of  the  control  system.  The  lower  frequency 
modes  have  larger  deflections  associated  with  them  and  are 
more  easily  controlled  than  the  higher  frequency  modes. 

Mode  24,  one  of  the  higher  frequency  modes,  indicated  that 
if  enough  sensor/actuator  pairs  are  affected  by  the  deflec¬ 
tion  associated  with  the  mode  shape,  active  damping  can  be 
successfully  employed.  Similarly,  in  the  lower  frequency 
modes,  the  nvimber  of  sensor/actuator  pairs  involved  affects 
the  damping  factor.  Modes  12  and  13  have  basically  the 
same  kind  of  mode  shape  (rotation  about  a  point  in  the 
upper  mirror  support  truss)  but  they  rotate  in  different 
planes.  Mode  12  has  one  sensor/actuator  pair  in  the  lower 
mirror  support  truss  to  control  the  vibration  and  had  a 
damping  factor  aroxind  .185.  Mode  13  affected  two  pairs 
in  the  lower  mirror  support  truss  and  had  a  damping  factor 
of  .22.  The  addition  of  more  control  pairs  in  the  lower 
mirror  support  truss  oriented  to  sense  motion  in  the  x  and 
y  directions  would  very  likely  increase  the  damping  of 
those  two  modes. 

It  seems  logical  that,  as  the  magnitude  of  a  pertur¬ 
bation  is  Increased,  the  damping  factor  would  show  a  con¬ 
sistent  trend,  either  up  or  down.  It  also  seems  likely  that, 
if  the  trend  in  the  change  of  the  damping  factor  is  not 
consistent,  there  should  be  some  point  where  the  damping 


factor  reaches  a  maxlmum/mininiuin  value  and  then  goes  down/ 
up  as  the  magnitude  of  the  perturbation  changes.  These 
trends  are  absent  in  much  of  the  data  developed  in  this 
research..  Possible  explanations  for  this  include  Inac- 
curacies  in  the  computer  modeling  of  the  structure  and 
the  extreme  flexibility  of  the  structure.  More  research 
on  the  structural  characteristics  of  the  CSDL  2  model  is 
required  if  this  model  is  to  be  used  to  develop  methods 
of  specifying  the  actual  control  system  hardware  for  a 
large  space  structure. 

Recommendations 

Several  areas  need  to  be  researched  in  detail  to  con¬ 
firm  the  results  of  this  thesis.  The  first  approach  that 
should  be  taken  is  to  determine  if  the  discrepancy  in  the 
modeling  of  the  two  stiffened  mirror  support  beams  has  an 
effect  on  the  results .  The  tubular  truss  should  be  modeled 
as  a  truss  and  not  as  an  equivalent  circular  beam.  This 
can  be  done  by  designating  appropriate  GRID  points  and 
"building"  a  beam  in  the  NASTRAN  input  deck.  Damping 
factors  in  this  remodeled  structure  should  then  be  com¬ 
pared  to  those  discussed  in  this  thesis. 

Secondly,  the  doubts  about  the  suitability  of  the  VCOSS 
version  of  this  model  should  be  resolved.  The  CSDL  2, 
Revision  1  model  should  be  analyzed  and  then  perturbed  in 
a  manner  similar  to  that  done  in  this  thesis.  This  will 
indicate  if  the  extreme  flexibility  in  the  VCOSS  version 
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of  this  model  (Revision  3)  is  responsible  for  the  erratic 
damping  behavior  discovered  in  this  research. 

Third,  the  effect  of  changing  the  orientation  of  the 
control  set  and  the  number  of  sensor/actuator  pairs  should 
be  studied,  with  emphasis  on  determining  a  control  system 
configuration  that  can  damp  the  modes  that  are  not  damped 
by  the  current  control  system. 

Finally,  the  predictability  of  the  effects  of  small 
perturbations  must  be  improved.  In  order  to  accomplish 
this,  the  data  base  must  be  expanded.  There  are  an  infi¬ 
nite  number  of  possible  perturbations;  this  research  has 
looked  at  46  of  them.  Several  of  these  perturbations 
indicate  that  small  deviations  in  mass  or  stiffness,  if 
they  exist  in  the  wrong  place,  can  negate  the  active  damp¬ 
ing  in  a  vibrational  mode  CFigures  33  and  36) .  Other 
cases  have  shown  that  the  damping  in  some  modes  is  in¬ 
creased  by  small  deviations.  The  predictability  of  these 
occurrences  must  be  improved  if  large  space  structures  are 
to  be  built  at  an  affordable  cost. 

Large  space  structures ,  whether  in  the  form  of  a 
telescope,  a  solar  power  station,  or  manned  space  stations, 
are  becoming  increasingly  important  in  the  plans  for  develop¬ 
ment  of  space.  In  order  to  be  able  to  afford  these  struc¬ 
tures,  the  specifications  for  their  components,  be  they 
control  systems  or  epoxy-graphite  tubes,  must  be  developed 
to  reflect  the  minimum  precision  required  to  accomplish 
the  mission.  This  thesis  has  explored  the  sensitivity  of 


a  LSS  to  small  deviations  in  mirror  mass,  beam  stiffness, 
and  mass  distribution.  Space  structures,  once  erected,  may 
exhibit  poor  controllability  due  to  small  manufacturing 
deviations  in  their  components  or  inaccuracies  in  the  design 
or  specification  of  the  control  system.  This  thesis  has 
shown  that  the  concept  of  "fine-tuning"  a  structure  in 
orbit  is  valid,  at  least  in  the  engineering  sense.  The 
ability  to  control  a  large  space  structure  is  dependent 
on  knowledge  of  the  interactions  between  active  damping, 
passive  damping,  bending,  twisting,  rotating,  and  the 
inertial  characteristics  of  the  large  masses.  Once  these 
are  understood,  the  large  space  structure  will  be  ready  to 
support  the  development  of  space. 
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Appendix  A.  NASTRAN  Bulk  Data  Listing 
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0.2018832-05 

100 
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100 
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PSAR» 
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6.1072562-04 

3.5617522-06* 

135 

*  135 

PBAR* 

3.5617522-06 

136 

7.1235042-06 

100 

6.1072562-04 

3.561  7522-06* 

136 

•  136 

P5AR* 

3.5617522-06 

137 

7.1235042-06 

100 

6.1072562-04 

3.5617522-06* 

137 

•  137 

PBAR* 

3.5617522-06 

138 

7.1235042-06 

100 

6.1072562-04 

3.5617522-06* 
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Appendix  B .  Computer  Methodology 

The  damping  ratio,  used  as  the  performance  index  in 
this  research,  is  the  end  result  of  a  lengthy  computational 
process.  The  programs  used  Include  NASTRAN,  FORMAT,  SELECT, 
ZETA,  and  D0FB34.  NASTRAM  is  widely  used  and  well  docimxented 
CRef  4)  and  will  not  he  discussed  here.  FORMAT,  SELECT,  ZETA, 
and  D0FB34  were  developed  at  the  Air  Force  Institute  of  Tech¬ 
nology  to  aid  in  the  analysis  of  the  CSDL  2  model.  Various 
forms  of  these  programs  exist;  the  users  must  ascertain  that 
they  use  the  appropriate  ones.  D0FB34  is  documented  and 
validated  in  reference  6.  These  programs,  as  listed  in 
the  following  appendices,  were  used  on  the  CDC  170/750  com¬ 
puter  system  at  Wright-Patterson  AFB,  OH.  The  first  task 
in  the  generation  of  data  was  to  run  NASTRAN  for  the  finite 
element  analysis  of  the  CSDL  2  model. 

NASTRAN  was  run  using  the  NASTJCL  control  cards  listed 
in  Appendix  C.  The  data  deck,  called  BULKl  in  NASTJCL,  con¬ 
tains  the  finite  element  model  of  the  LSS.  It  is  listed  in 


Appendix  A.  The  output  from  the  NASTRAN  program  is  filed 
as  ACOSSFUN  and  contains  a  six  line  NASTRAN  header,  followed 
hy  the  eigenvalue,  mode  number,  and  matrix  for  each  mode. 
The  entry  in  line  three  of  the  header,  titled  LABEL,  is  used 
to  keep  track  of  the  particular  bulk  data  deck  used  as  input 
to  NASTRAN.  This  variable,  labeled  DATASET  in  succeeding 


programs,  is  passed  through  the  rest  of  the  computer  analysis 
and  printed  on  the  final  product. 

Program  FORMAT  (Appendix  D)  is  used  to  reformat  ACOSSPUN 
into  the  format  required  for  the  SELECT  program.  FORMAT  reads 
ACOSSPUN  C^elabeled  TAPES)  and  prints  out  (on  TAPES)  the  nxim- 
ber  of  critical  modes  C12  in  this  case)  and  the  specific  mode 
numbers  to  be  analyzed  (.4,5 .6, 7, 12 , . . .  29) ,  followed  by  DATASET, 
the  mode  number,  eigenvalue,  and  (p  matrix  for  each  mode.  This 
program  should  be  compiled  using  the  FORTRAN  IV  compiler. 

SELECT  (Appendix  E)  reads  TAPES  and  then  outputs  DATASET, 
the  matrix  for  each  of  the  12  modes,  and  the  12  eigenvalues 
onto  TAPES.  Next,  the  number  of  modes  being  analyzed,  the 
number  of  actuators,  the  nvimber  of  sensors,  and  the  location 
and  orientation  of  the  sensors/actuators  are  added  to  the 
front  of  TAPES.  An  example  of  this  set  of  data  is  listed 
in  Appendix  G. 

Program  ZETA  (Appendix  F)  prepares  the  data  for  input 
to  D0FB34.  ZETA  reads  the  revised  TAPES  and  computes  the  B 
and  C  matrices.  It  also  takes  the  square  root  of  the  eigen¬ 
value  to  get  the  natural  frequency  in  radians  per  second. 

As  discussed  earlier,  the  B  and  C  matrices  are  the  transpose 
of  each  other  if  the  sensors  and  actuators  are  collocated. 
DOFB34,  the  next  program  used  in  the  analysis,  requires  C 
in  its  transposed  form  and  then  transposes  it  back  to  its 
original  form.  The  B  matrix  is  therefore  printed  by  ZETA 


to  represent  the  transposed  C  matrix.  Various  other  param¬ 
eters  for  D0FB34  are  also  printed  by  ZETA  on  an  output  file . 
These  parameters  are  explained  In  reference  6 .  Thus ,  the 
output  from  ZETA  contains  DATASET,  2  lines  of  parameters  for 
DOFB34,  the  B  matrices  for  the  12  modes,  the  C  matrix  for 
the  12  modes  Crepresented  by  B  In  this  case) ,  the  12  natural 
frequencies,  and  finally,  9  lines  of  parameters  for  D0FB34. 

DOFB34  reads  the  output  file  from  ZETA,  relabeled  TAPES, 
and  calculates  the  complex  conjugate  pairs  that  are  used  to 
calculate  the  damping  factor  for  the  various  modes.  These 
pairs  are  titled  "Eigenvalues  Al+CBSTARl) (KSTARl) (Cl) ,  A2+ 
(BSTAR2) CKSTAR2) (C2) ,  and  A3+CBSTAR3) CKSTAR3) (C3) . "  Each 
of  the  12  modes  Is  represented  by  two  eigenvalue  pairs  having 
the  same  real  value  and  only  different  signs  on  the  imaginary 
value.  Thus,  the  original  NASTRAN  data  has  been  reduced  to 
the  24  eigenvalue  pairs  output  by  D0FB34. 

The  damping  factor  calculation  (equation  14)  uses  the 
absolute  value  of  the  real  and  imaginary  values. 
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Appendix  C .  NASTRAN  Job  Control  Listing  CNASTJCL) 


The  commands  listed  below  are  used  to  run  the  NASTRAN 
program.  The  job  control  card  should  include  a  time  limit 
of  at  least  8QQ  seconds,  an  Input/output  limit  of  1000 
seconds,  and  170Q0Q  xmits  of  core  memory.  The  NASTRAN  input 
deck  is  labeled  BULKl  in  this  listing,  and  the  NASTRAN  output 
(pxjnch  file)  goes  to  the  file  called  PUN.  ACOSSPUN  is  the 
permanent  file  name  assigned  to  PUN  during  filing . 


C.Job  Control  Card) 

ATTACH. NASTRAN ,  NASTRAN . ID-NASTRAN , SN-AFFDL ,MR-1 . 
ATTACH . NASI . NASI , ID-TRAN . SN-ASDAD . 

ATTACH , NAS2 . NAS2 , ID-TRAN , SN-ASDAD . 

LIBRARY, NAS 1,NAS2. 

ATTACH, DATA, BULKl. 

REQUEST, PUN, *PF. 

RFL, 17000. 

NASTRAN, DATA, PUN.  ATTACH 
REWIND, PUN. 

CATALOG , PUN , ACOSSPUN , CY-1 , RP-999 . 

★EOF 
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Appendix  D.  Program  FORMAT 
The  first  two  WRITE  statements  write  the  niunber  of 
modes  to  be  selected  and  the  particular  modes  that  make  up 
the  selection  onto  the  output,  TAPES.  This  data  is  required 
by  SELECT.  The  program  then  searches  for  the  name  "TITLE" 
in  the  NASTRAN  output  file  Ccalled  ACOSSPUN  and  relabeled 
TAPES)  by  assigning  the  variable  ANAME  to  the  first  ten  let¬ 
ters  of  every  line  and  continuing  when  ANAME>TITLE.  It  then 
reads  SUBTI,  which  is  the  subtitle,  and  DATASET,  the  variable 
assigned  to  the  LABEL  block  in  the  NASTRAN  input  deck.  Af¬ 
ter  reading  the  next  three  lines  (not  used  for  anything  here) , 
it  reads  EIGEN  Cthe  eigenvalue  of  the  mode)  and  MODEN  (the 
mode  number)  and  writes  DATASET,  MODEN,  and  EIGEN  on  TAPES. 

The  matrix  for  each  of  the  modes  CT1,T2,T3,R1,R2,R3)  is 
then  read  and  written  onto  TAPES.  This  process  continues 
until  the  data  for  all  of  the  modes  has  been  reformatted 
and  written  onto  TAPES. 


PROGRAM  FORMAT (INPUT , OUTPUT , TAPES , TAPES ) 
DIMENSION  CARDCS) 

WRITE (S, 11)  12 

WRITE(8,12)  4, S. 6, 7, 12, 13, 17, 21,22,24,28,29 

11  F0RMATCI2) 

12  FORMAT (12 (12, IX)) 

1  READ (S,  100)  ANAME 

IF (EOF (S)) 1000, IS 
100  FORMAT (AlO) 

IS  IF  (ANAME. EQ.10H$TITLE  -)  GO  TO  ISO 
GO  TO  1 

ISO  READ(5,100)  SUBTI 
READ(S,23S)  DATASET 
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DO  200  1-1,3 

READC5 , 205)  (CARDCJ)  .  J»1 , 8) 
205  FORMATC8A10) 

200  CONTINUE 
10  FORMAT (13, E13. 6) 

235  FORMAT (:10X,2A10) 

READC5,210)  EIGEN, MODEN 
210  F0RMATC14X,E14.6.11X.I3> 
WRITE C8, 235)  DATASET 
WRITE (8, 10)  MODEN, EIGEN 
DO  220  1-1,59 
READC5,215)  T1,T2,T3 
215  F0RMATC23X,3(E13.6,5X)) 

•  READC5,215)  R1,R2,R3  ' 

WRITEC8,225)  T1,T2,T3 
225  FORMAT C3E15. 6) 

WRITEC8,225)  R1,R2,R3 
220  CONTINUE 
GO  TO  1 
1000  STOP 
END 
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Appendix  E.  Program  SELECT 

This  program  reads  the  TAPES  prepared  by  FOKMAT.  The 
variables  are  defined  as  follows : 

N  is  the  number  of  modes  to  be  selected. 

ID  (I)  is  a  matrix  of  the  mode  nvimbers  to  be  selected. 

DATASET (K)  is  a  matrix  of  the  bulk  data  identifier. 

MODECK)  is  the  matrix  of  mode  numbers. 

EIGEN (K)  is  the  matrix  of  eigenvalues. 

MAT(I,J,K)  is  a  59  by  6  X  30  matrix  of  the  41  matrices. 
DATASET  is  printed  onto  TAPE6,  followed  by  the  4)  matrix  for 
the  selected  modes.  Finally,  the  eigenvalues  of  the  selected 
modes  are  listed.  The  MODES  READ  and  MODES  FILED  statements 
are  displayed  on  the  output  device  (the  screen)  to  assure  the 
user  that  the  proper  actions  have  been  accomplished. 


PROGRAM  SELECT (INPUT. OUTPUT, TAPES, TAPE6) 

DIMENSION  MAT(59,6.30) .MODE(30) ,ID(30) ,DATASET(30) 
DIMENSION  EIGEN(30) 

CHARACTER*40 , DATASET 
READ(8.*)N 

READ(8.*)(ID(I) .1-1. N) 

DO  1  K-1,30 

READ (8 . ’ (A40) * )  DATASET (K) 

READ(8,100)MODE(K) ,EIGEN(K) 

READ(8,200) ((MAT(I,J,K) ,J-1, 6) .1-1,59) 

1  CONTINUE 

WRITE (6 , ’ (A40 ) ' )  DATASET ( 1) 

DO  2  K-1,N 
M-ID(K) 

WRITE (6, 300) ((MAT (I, J.M) ,J-1, 6) ,1-1,59) 

2  CONTINUE 
DO  3  K-1,N 
M-ID(K) 

WRITE(6,400)  EIGEN(M) 

3  CONTINUE 
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PRINT*, 'MODES  READ  ARE  ; 
PRINT*, 'MODES  FILED  ARE 
100  FORMAT (13, E13. 6) 

200  FORMAT (3E15. 6) 

300  FORMAT (3E 15. 6) 

400  FORMAT (E13. 6) 

END 


• , (MODE (I) ,1-1,30) 
',(ID(I),I-1,N) 


I 


This  program  reads  the  TAPE6  output  from  SELECT  with  the 
control  system  orientation  added  onto  the  front  of  TAPES. 

The  variables  are  defined  as  follows: 

N  is  the  nxmiber  of  modes  to  be  analyzed. 

NA  is  the  number  of  actuators. 

NS  is  the  nimber  of  sensors. 

ALPHACI) ,  BETA(I) ,  and  GAMMA(I)  are  the  orientation  of 
the  actuators. 

NODE (I)  is  the  location  (node  point)  of  the  actuators. 
ALPHS(I),  BETS(I).  GAMMS(I)  and  NODS(I)  represent  the 
orientation  and  location  of  the  sensors . 

ID (I)  is  a  59  unit  array  of  the  node  point  numbers  of 
CSDL  2.  (This  is  added  onto  TAPES  with  the  control 
system  data.) 

PHI (354, 12)  are  the  (j)  matrices  of  the  12  selected  modes. 
(Note  that  59  nodes  times  S  values  per  node  equals  354 
values . ) 

B(12,21)  and  C (2 1,12)  are  the  B  and  C  matrices  calculated 
in  this  program. 

FREQ(I)  is  the  natural  frequency  in  radians  per  second, 
it  is  the  square  root  of  EIGEN(I) , 

EIGEN (T)  and  DATASET  remain  the  same  from  SELECT. 

After  reading  in  the  revised  TAPES,  the  six  angles  (ALPHA, 
BETA,  etc.)  are  converted  from  degrees  to  radians  by  CON. 


Next,  the  program  matches  up  the  node  numbers  and  the 
actuator/sensor  locations,  finds  the  proper  part  of  the  41 
matrix  using  K  *  6*K-5 ,  and  calculates  the  B  and  C  matrices 
using  equation  (13)  for  B  and  the  transpose  of  that  for  C. 
After  writing  DATASET,  2  lines  of  input  for  D0FB34  are 
written,  followed  by  the  B  matrix  twice  (D0FB34  wants  and 
■  B) .  This  output  is  not  automatically  written  onto  a 
TAPE;  the  disposition  of  the  output  must  be  specified  when 
the  program  is  rxm.  After  taking  the  square  root  of  EIGEN, 
FREQ  is  written.  Finally,  the  other  parameters  for  D0FB34 
are  written  and  the  program  ends . 


PROGRAM  ZETA (INPUT , OUTPUT , TAPES , TAPES ) 

DIMENSION  PHI (354, 12) ,ALPHA(21) ,BETA(21) , GAMMA(21) 

DIMENSION  ALPHS(J21)  ,BETS(.21)  , GAMMS (21)  ,N0DE(21)  ,N0DS(21) 

DIMENSION  ID(59) ,B(12,21) ,C(21,12) , EIGEN ( 12) . FREQ (12) 

CHARACTER*40 , DATASET 

READ(6,*)N,NA,NS 

READ (6,*) (ALPHA (I) ,1-1, NA) 

READ (6,*) (BETA(I),I-1,NA) 

READ(6,*) (GAMMA(I) ,I-1,NA) 

READ (6,*)  (NODE (I)  ,1-1, NA) 

READ(6,*)(ALPHS(I) ,1-1, NS) 

READ (6,*)  (BETS (I)  ,1-1, NS) 

READ(6,*) (GAMMS(I) ,1-1, NS) 

READ (6,*) (NODS (I) ,1-1, NS) 

READ(6,*) (ID(I) ,1-1,59) 

READ (6 . * ' (A40) ’ )  DATASET 

READ (6, 400) ((PHI (I, J) ,1-1,354) ,J-1,N) 

READ(6,230) (EIGEN (K) ,K-1,N) 

PI-3.14159265E0 
CON-PI/ 18 0.0 
DO  30  J-1,NA 
ALPHA (J) -ALPHA (J) *C0N 
BETA ( J) -BETA ( J ) *CON 
GAMMA(J) -GAMMA (J) *CON 
DO  10  K-1,59 

10  IFCID(K) .EQ.N0DE(J))G0  TO  20 

20  CONTINUE 
K-6*K-5 
DO  30  r-l,N 
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B (I , J) =PHI (K . I) *COS (ALPHA( J) ) +PHI (K+1 , I) *COS (BETA( J) ) 

1  +PHI(K+2.1)*C0S(GAMMA(J)) 

30  CONTINUE 
DO  60  1=1, NS 
ALPHS ( I ) =ALPHS ( I ) *CON 
BETS(I)=BETS(I)*CON 
GAMMS ( I ) -GAMMS ( I) *CON 
DO  40  K-1,59 

40  IF  (ID(K) .EQ.NODS(I))GO  TO  50 
50  CONTINUE 
K-6*K-5 
DO  60  J-1,N 

C (I , J) -PHI (K , J) *COS (ALPHS (I ) ) +PHI (K+1 , J) *COS (BETS (I ) ) 

1  +PHI(K+2.J)*COS(GAMMS(I)) 

60  CONTINUE 

WRITE (* , ' (A40) ' )  DATASET 
PRINT' ("3") ' 

PRINT' ("3  4  5  0  21  21  0.01")’ 

DO  500  1=1. N 

500  PRINT  220, (B(I,J) ,J=1.NA) 

DO  600  1=1, N 

600  PRINT  220, (B(I.J) .J=1,NS) 

DO  700  I-1,N 

FREQ ( I ) -SQRT (EIGEN ( I ) ) 

PRINT  230.FREQ(I) 

700  CONTINUE 

PRINT' ("1  2  3")’ 

PRINT' ("4  5  6  7")' 

PRINT' ("8  9  10  11  12")' 

PRINT' ("0") ' 

PRINT' ("0") ' 

PRINT' ("2  2  2  5000  5000  5000  5000  5000  5000  5000  5000  5000")' 
PRINT' ("2") ' 

PRINT' ("2") ' 

PRINT' ("3") ' 

220  FORMAT(lX,/,6(/.6X.4Fl0.6)) 

230  FORMAT(E13,6) 

400  FORMAT (3E15. 8) 

END 
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Appendix  H.  B  Matrix  for  the  Unperturbed  Model 


This  appendix  shows  the  B  matrix  conqjuted  by  program 


ZETA  for  the  unperturbed  model.  B  is  a  12  by  21  matrix; 
the  printout  lists  the  row  number,  followed  by  the  21 
elements  in  that  row. 
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Appendix  I .  Mode  Shapes  of  the  Unperturbed  Model 

Figure  15  through  24  show  the  mode  shapes  for  the  12 
modes  under  consideration  in  this  research.  These  plots 
were  made  using  the  GCSNAST  plotting  program. 
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Figure  18.  Mode  Shape  13 
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Figure  23.  Mode  Shape  28 


Appendix  J .  Damping  Factors 


This  appendix  contains  the  actual  damping  factors  for 
the  12  modes  and  the  various  perturbations.  Modes  4,  5, 
and  6  alvays  had  a  damping  factor  of  zero.  The  damping 
factor  for  mode  7  was  always  greater  than  one.  Modes 
22,  28,  and  29  had  damping  factors  of  .01000  throughout 
the  research.  The  values  for  the  remaining  modes  are  listed 
below. 


Table  V.  Damping  Factors  for  the  Unperturbed  Model 
Mode  Damping  Factor 


12 

.04386 

13 

.22070 

17 

.18471 

21 

.01013 

24 

.01159 
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Table  VI.  Damping  Factors  for  Mirror  Mass  Changes 


(Effected 

Mirrors 

% 

Change 

Damping  Factor  For  Mode 

12 

13 

17 

21 

24 

All 

-10 

.18910 

.19923 

.04902 

.01016 

.01188 

-  5 

.17009 

.21432 

.04771 

.01020 

.01231 

-  1 

.15497 

.21392 

.04877 

.01015 

.01176 

+  1 

.18453 

.22180 

.04392 

.01015 

.01183 

+  5 

.18355 

.22611 

.04413 

.01015 

.01178 

+10 

.17654 

.21962 

.04850 

.01013 

.01154 

2  and  3 

-10 

.16899 

.20711 

.05167 

.01015 

.01184 

-  5 

.17130 

.21043 

.05210 

.01018 

.01164 

-  1 

.17160 

.21691 

.04796 

.01018 

.01210 

+  1 

.17259 

.21759 

.04808 

.01022 

.01254 

+  5 

.17462 

.21851 

.04831 

.01016 

.01183 

+10 

.18685 

.22607 

.04434 

.01020 

.01237 

2 

-10 

.16720 

.21154 

.04753 

.01021 

.01247 

-  5 

.19068 

.20516 

.04935 

.01020 

.01228 

-  1 

.18428 

.22024 

.04380 

.01019 

.01120 

+  1 

.18516 

.22109 

.04391 

.01015 

.01174 

+  5 

.19900 

.20527 

.05006 

.01017 

.01198 

+10 

.18900 

.22381 

.04437 

.01015 

.01179 

3 

-10 

.17394 

.21279 

.05239 

.01014 

.01174 

-  5 

.19421 

.20608 

.04974 

.01017 

.01199 

-  1 

.18485 

.22053 

.04386 

.01014 

.01170 

+  1 

.17207 

.21724 

.04804 

.01017 

.01195 

+  5 

. 19446 

.20604 

.04969 

.01016 

.01183 

+10 

.19451 

.20605 

.04965 

.01022 

.01252 

114 


Table  VII.  Damping  Factors  for  Stiffness  Increases 
(See  Figures  4-12) 


Perturbation 

%  Change 

Damping 

Factor  for  Mode 

Set 

in  Area 

12 

13 

17 

21 

24 

1 

+10 

.19415 

.20606 

.04967 

.01018 

.01217 

+50 

.17090 

.21740 

.04871 

.01015 

.01215 

2 

+10 

.18406 

.22113 

.04382 

.01014 

.01174 

+50 

.19358 

.20590 

.04950 

.01013 

.01194 

3 

+10 

.17376 

.21252 

.05249 

.01015 

.01174 

+50 

.19461 

.20598 

.04970 

.01015 

.01182 

4 

+10 

.17375 

.20471 

.05249 

.01013 

.01161 

+50 

. 19446 

.20605 

.04970 

.01016 

.01192 

5 

+10 

.19435 

.20606 

.04972 

.01016 

.01182 

+50 

.19436 

.20606 

.04972 

.01016 

.01187 

6 

+10 

.19435 

.20606 

.04972 

.01016 

.01192 

+50 

.19437 

.20605 

.04972 

.01019 

.01223 

7 

+10 

.16815 

.21613 

.01004 

.01018 

.01215 

+50 

.17422 

.21292 

.05254 

.01018 

.01213 

8 

+10 

.16972 

.22219 

.01070 

.01014 

.01173 

+50 

.17238 

.21774 

.04804 

.01020 

.01256 

9 

+10 

.17293 

.21729 

.04802 

.01020 

.01231 

+50 

.17200 

.21733 

.04801 

.01019 

.01234 

Table  VIII.  Damping  Factors  for  Midsection  Mass  Increases 


(See  Figure  13) 


Increase  in 
nass  per  bar,  kg 

Damping 

Factor  for  Mode 

12 

13 

17 

21 

24 

5 

.17216 

.21734 

.04801 

.01016 

.01195 

25 

.16889 

.21527 

.01003 

.01018 

.01222 

50 

.18547 

.22218 

.04378 

.01018 

.01215 

85.9 

.20089 

.20314 

.04951 

.01019 

.01226 

Appendix  K.  GCSNAST  Procedures 


The  GSCNAST  plotting  program  was  used  to  produce  the 
deflected  view  of  the  CSDL  2  model  in  Appendix  I.  GCSNAST 
uses  the  NASTRAN  hulk  data  deck  for  data  on  the  shape  and 
size  of  the  structure  and  then  refers  to  a  deflection  file 
that  is  processed  from  the  NASTRAN  output  for  data  on  the 
shape  of  the  deflected  structure.  The  first  step  in  running 
GCSNAST  is  to  create  the  deflection  file.  The  following 
commands  accomplish  this: 

ATTACH , NASTPP , NASTPP , ID-GCSNAST , SN-AFFDL . 

ATTACH , NAS 1 . NAS 1 , ID-TRAN . SN-ASDAD . 

ATTACH . NAS2 , NAS2 , ID-TRAN , SN-ASDAD . 

LIBRARY, NASI, NAS2. 

ATTACH, PUN , ACOSSPUN , CY-1 

REQUEST , DEFORM , *PF . 

NASTPP, PUN. 

The  command  NASTPP, PUN  runs  the  program  and  puts  the  deflec¬ 
tion  file  into  the  file  labeled  DEFORM.  ACOSSPUN  is  the 
permanent  file  name  used  for  the  NASTRAN  output  (punch)  file. 

After  creating  the  deflection  file,  GCSNAST  can  be  run. 
The  following  two  commands  retrieve  the  GCSNAST  program  and 
set  up  its  files : 

ATTACH , GCSNAST , GCSNAST , ID-GCSNAST , SN-ASDAD . 

GCSNAST. 

Next,  the  program  asks  for  the  number  of  grid  points 
and  elements  in  the  model,  the  type  of  data  (NAST) ,  and, 
finally,  the  name  and  location  of  the  file  that  contains 
the  NASTRAN  input  deck.  At  this  point,  the  program  is  ready 
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to  use  and  the  user  can  refer  to  the  GCSNAST  manual  for  the 
specific  commands  used  in  GCSNAST.  The  figures  in  Appendix  I 
were  drawn  using  the  DEFORM  command  and  the  deflection  file. 
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this  large  space  telescope  were  varied  by  up  to  10%,  selected  structural 
elements  were  stiffened  by  increasing  their  cross-sectional  areas  by  10  and 
5.0%,  and,  finally,  two  structural  elements  in  the  middle  of  the  telescope 
were  stiffened  to  represent  the  addition  of  a  lumped  mass  located  away  from 
the  control  system  sensors  and  actuators.  A  control  system  of  21  collocated 
sensors  and  actuators,  positioned  at  the  top  and  bottom  of  the  telescope, 
was  used  in  this  analysis. 

The  analysis  was  accomplished  using  NASTRAN  for  the  finite  element 
analysis,  and,  after  selecting  certain  vibration  modes  for  further  study, 
the  complex  conjugate  pairs  used  to  determine  the  damping  factors  were 
calculated.  This  last  step  was  accomplished  through  spillover  reduction 
and  reduction  of  the  closed  loop  matrix  to  a  lower  triangular  form.  This 
procedure  was  repeated  for  each  of  the  46  perturbations  applied  to  the  model 
in  this  research.  The  results  indicate  that  the  damping  factors  for  the  CSDL 
model  are  quite  sensitive  to  small  perturbations  and  that  it  is  very  diffi¬ 
cult  to  predict  the  effect  that  a  perturbation  to  the  model  will  have  on  the 
damping  factors.  If  the  vibrations  in  an  LSS  are  to  be  altered  by  changing 
mass  or  stiffness,  a  thorough  analysis  is  required  to  ensure  that  increasing 
the  damping  for  one  mode  does  not  result  in  an  undesirable  change  in  the 
damping  for  another  mode. 
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